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We present predictions for tiie flux averaged muon energy spectra of quasi-elastic (QE) and 1-pion 
production events for the K2K long-baseline experiment. Using the general kinematical consid- 
erations we show that the muon energy spectra closely follow the neutrino energy spectrum with 
downward shift of the energy scale by 0.15 GeV (QE) and 0.4 GeV (1-pion production). These 
predictions seem to agree with the observed muon energy spectra in the K2K nearby detector. We 
also show the spectral distortion of these muon energy spectra due to the neutrino oscillation for 
the SK detector. Comparison of the predicted spectral distortions with the observed muon spectra 
of the 1-Ring and 2-Ring muon events in the SK detector will help to determine the oscillation 
parameters. The results will be applicable to other LBL experiments as well. 
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I. INTRODUCTION 

Recently the KEK to Kamioka long-baseline neutrino 
oscillation experiment (K2K) has published its first re- 
sult Q, which confirms the existence of t'^ oscillation 
as seen in the Super-Kamiokande (SK) atmospheric neu- 
trino data 0. The observed oscillation parameters from 
K2K agree well with the neutrino mass and mixing an- 
gles deduced from the atmospheric neutrino oscillation 
data 

sin^ 29 ~ 1 and Am' ~ 3 x IQ-^eV^ . 

As is well known, in a two flavor scenario, the prob- 
ability for a muon neutrino with energy to remain a 
muon neutrino after propagating the distance L is given 
by the following expression 
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Basically, the standard approach to measure the oscil- 
lation parameters is to determine the oscillation proba- 
bility in Eq. in dependence of E^. At the position of 
the minimum Am' can be determined from the condition 
= f and sin^ 29 from P^f,{E,^nun) = 1 - sin^ 29. 
The neutrino energy is not directly measurable but can 
be reconstructed from the simple kinematics of quasi- 
elastic (QE) scattering events. Measuring the energy E'^ 
and the polar angle cos 0^ of the produced muon allows to 
reconstruct E^, with help of the following relation (even 
if the scattered proton is not observed) 



= £;^[£:^, cos6'^] = 
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M - Ef, + \kf,\cost 
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Here M denotes the proton mass, the muon mass and 
is the three-momentum of the muon in the laboratory 
system. 



However, in practice there are some difficulties. First 
of all, the experimental one-ring muon events (IR/^) are 
not pure QE event samples. About 30% of the IR/x 
events are 1-pion production events with unidentified or 
absorbed pions. For the 1-pion events Eq. Q would sys- 
tematically underestimate the true neutrino energy j^]- 
Secondly, the reconstruction of E^, gets more complicated 
including binding energy es and Fermi motion of the tar- 
get nucleons 
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where p is the three momentum and Ep ~ ^ M' + p' 
the energy of the initial nucleon. Further, 9p is the polar 
angle of the target nucleon w.r.t. the direction of the 
incoming neutrino. Neglecting es and the momentum p 
Eq. is recovered. Since the momentum p is unknown, 
"Si \v\ < Pf where pf is the Fermi momentum, this 
will lead to an uncertainty of the reconstructed neutrino 
energy at given values E^, cos 9^, and es of about -9% 
to +6% for a single event. 

Hence we can see no reliable way for reconstructing 
the neutrino energy for the IR/i sample on an event by 
event basis. On the other hand the muon energy is a 
directly measurable quantity for each event. Therefore it 
seems to us to be a better variable for testing the spectral 
distortion phenomenon compared to the reconstructed 
neutrino energy. 

In this ta,lk we summarize the basic ideas and the main 
results in 4] where we have used kinematic considerations 
to predict the muon energy spectra of the QE and 1-pion 
resonance production events which constitute the bulk 
of the charged-current ly^^ scattering events in the K2K 
experiment. These predictions can be checked with the 
observed muon energy spectra from the nearby detector. 
We also present the distortion of these muon spectra due 
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to Vfj^ oscillation, which one expects to see at the SK 
detector. Comparison of the predicted muon spectra with 
those of the observed QE and 1-pion events at the SK 
detector will be very useful in determining the oscillation 
parameters. 



II. FLUX AVERAGED MUON ENERGY 
SPECTRA 

The flux averaged muon energy spectra for QE and 
1-pion events are given by 



(— > 



(4) 



where f{E^) is the neutrino flux at K2K for the nearby 
detector(ND) and 'R' denotes QE and A resonance con- 
tribution to 1-pion production, which dominates the lat- 
ter. Simple kinematic considerations lead to the following 
approximation for the flux averaged muon energy spec- 
tra, both, for QE and 1-pion production |Jj 
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with 



En = E^ + AE^ = E^ 



0.15 GeV for QE 
0.4 GeV for 1-pion 



(6) 



Furthermore, it is well known that the total cross sec- 
tions for QE and A production tend to constant values 
for neutrino energies of about 1 GeV and 1.4 GeV, re- 
spectively: a^f[E,y] — » A''-^. Hence, for muon energies 
larger than about 1.2 GeV, Eq. can be further sim- 
plified by replacing ct/J^ by its constant asymptotic value 



for Ef, > 1.2 GcV , (7) 



with 



_ j 4.5 fb for QE 



5.5 fb for 1-pion 



(8) 



The normalizations correspond to the average cross- 
section per nucleon for a H2O target 0]. 

Thus we conclude that at large muon energies E^ > 
1 GeV the flux averaged muon energy cross section {^-) 
is directly proportional to the neutrino flux shifted in en- 
ergy. The normalizations N^- and energy shifts AE^ are 
predictions of our theoretical calculation Q which can be 
verified experimentally with the muon energy spectra of 
the QE and 1-pion events observed at the ND. Further- 
more, Eqs. lO, and also apply to the far detector 
(FD) if one replaces the flux at the ND by the flux at the 



FD which is distorted by the neutrino oscillation proba- 
bility Pf,f,{E,): 



f{Eu) ^ kK{E,^ = f{E,) X P^^iE,) 



(9) 



Comparing these predictions with the observed muon en- 
ergy spectra of the QE and 1-Pion events of the SK detec- 
tor will test the spectral distortion due to i'^ oscillation 
and determine the oscillation parameters. Particularly, 
on the higher energy side of the peak the relative size of 
the SK to the ND cross-sections provides a direct mea- 
sure of the spectral distortion and hence the underlying 
oscillation parameters: 
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In Sec, mil we present exact calculations of the QE [g 
and 1-pion production cross sections. The dominant con- 
tribution to 1-pion production is from A resonance pro- 
duction for which we take the formalism of Refs. 7, 8, Qj. 
For completeness we have also included the Pii(1440) 
and 5*11(1535) resonance contributions for which we used 
the parameterization and the form factors of Ref. [lol| . 
We estimate the contribution from the still higher reso- 
nances along with the non-resonant background to be no 
more than 5 — 10% of the 1-pion production cross section 
at these low energies (see Fig.^. Therefore the accuracy 
of our prediction should be as good as that of the K2K 
experiment. 

We compare our exact calculations with the approx- 
imation in Eqs. (0 and H10|l from which the range of 
validity can be inferred. A comparison with the approx- 
imation in Eq. Q can be found in Ref. Q (see Figs. 2 
and 3). 



III. NUMERICAL RESULTS 

In this section we present our numerical results. Fig. 
n shows the K2K neutrino energy spectrum. In addi- 
tion, we show an approximated spectrum by a Lorentzian 
given by 



ME. 
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TT {E,-Eo)^+-f^ 
Eo = 1.2 GeV, r = 0.6 GeV 



(11) 



where TV is an appropriate normalization factor. 

Fig. shows the predicted muon energy spectra for 
the QE (solid line) and 1-pion production (dashed line) 
processes. Clearly one can see that the peak at Eq = 
1.2 GeV is shifted left by AE ~ 0.15 GeV for the QE 
and AE ~ 0.4 — 0.5 GeV for the A resonance produc- 
tion which dominates the 1-pion production process. The 
steepness of the muon energy spectra at the low energy 
reflects the threshold rise of and the steep neutrino 
flux. On the other hand the muon energy spectra closely 
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FIG. 1: K2K neutrino energy spectrum. The solid line is the 
exact spectrum normalized to unit area. The dashed line shows 
the approximated spectrum by the Lorentzian in Eq. Illlll . We 
also show for comparison as dotted line the same Lorentzian, but 
with normalization 1.4 instead of 1.25. 



follow the shape of the neutrino energy spectrum on the 
right side of the peak. The predicted exact muon energy 
spectra of Fig. agree reasonably well with the corre- 
sponding spectra of the K2K ND ij for both the QE and 
the non-QE sample. In particular one can compare the 
predicted QE spectrum with their simulated QE spec- 
trum shown in Fig. 1 of Ref. 0. Their Figs, la and c 
show separately the QE muon momentum distribution 
for the 1-Ring muon (IR/i) sample of the IKT and the 
QE enhanced sample of the FGD respectively. The two 
play complementary roles in covering the complete muon 
energy range, as the IKT and the FGD have high effi- 
ciencies at Efj^ < 1 GcV and E^j, > 1 GeV, respectively 
0. One can not compare our predicted muon energy 
spectra with these figures quantitatively without folding 
in these efficiency factors, which are not available to us. 
But there is good qualitative agreement between the pre- 
dicted QE spectrum of our Fig. with their Fig. Ic at 
Ef, > 1 GeV and Fig. la at < 1 GeV. While the 
former shows the position of the peak and the shape of 
the spectrum to the right, the latter shows the broaden- 
ing of spectrum down to ~ 0.4 GeV. Similarly one 
sees good agreement between the predicted muon energy 
spectrum of our Fig. [2^ for 1-pion events with the non- 
QE spectra of their Fig. lc,d at E'^ > 1 GeV and Fig. 
la at Ef^ < 1 GeV. Thus one has a simple and robust 
prediction for the shape of the muon energy spectrum 
in terms of the neutrino spectrum not only for the QE 
events but also for the 1-pion production events, which 
dominate the inelastic events. 

Fig. I^b shows the corresponding muon energy spectra 
of the QE and 1-pion events for the SK detector, pre- 
dicted by Eqs. Q and ©• One can clearly see the dis- 
tortion of the muon energy spectrum due the oscilla- 
tion. They should be compared with the observed muon 
energy spectra of the 1-Ring and 2-Ring muon events at 
the SK detector, after taking into account the pion detec- 
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FIG. 2: Exact predictions of the muon energy spectra for the (a) 
Nearby and (b) SK detectors of the K2K experiment. The QE 
(solid line) and the 1-pion production (dashed line) cross-sections 
are shown along with the A resonance contribution (dotted line). 



tion efficiency. We hope such a comparison will be done 
by the K2K collaboration. 

For the estimation of the pion detection efficiency it 
is necessary to consider nuclear effects, i.e., Pauli block- 
ing, nuclear absorption and charge exchange of the pro- 
duced pions which can rescatter several times in the nu- 
cleus. Therefore we have included these effects following 
the prescription of Refs. |5j|, Jl^ J^J. Since the dominant 
contribution to 1-pion production processes comes from 
A resonance production on oxygen, we have evaluated 
the effects of nuclear absorption and rescattering on the 
produced pions for this case. The relevant charged cur- 
rent subprocesses are vp — > fj.^pir'^, vn — > /i^nvr^ and 
vn fi^pn'^ with relative cross-sections 9:1:2 for the 
dominant contribution from A resonance. 

Fig. shows the effects of nuclear corrections on the 
produced 7r+ and 7r° spectra from these processes for 
the nearby detector averaged over the neutrino spectrum. 
The results are very similar for the SK detector. Nuclear 
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rescattering effects result in enhancing the tt" events at 
the cost of the dominant 7r+ component. But taken to- 
gether we see a nearly 20% drop in the rate of 1-pion 
events due to nuclear absorption of the produced pion. 
Moreover about 10% of the remaining events corresponds 
to the pion momentum being less than the Cerenkov 
threshold of 100 MeV. Therefore one expects about 70% 
of the A events to give a detectable pion ring at the 
SK detector while the remaining 30% appears like a QE 
event. Adding the latter to the 35% of genuine QE events 
would imply that about 50% of the CC events will appear 
QE-like at the SK detector. Alternatively the observed 
muon energy spectrum of the sum of 1-Ring and 2-Ring 
muon events could be compared with the predicted spec- 
trum of the sum of QE and 1-Pion events. Although a 
part of the 2-Ring events may come from multi-pion pro- 
duction, the resulting error may be small since multi-pion 
events at the ND constitute only ^ 15% of CC events. 
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FIG. 4: The predicted exact muon energy spectra (solid lines) 
and the approximation (dashed lines) according to Eq. ||7|l for 
QE and 1-pion production for the (a) Nearby detector and (b) 
SK detector of the K2K experiment. 
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FIG. 3: The momentum distribution of the decay pion for 
charged current resonance production in oxygen with (solid line) 
and without nuclear correction (dotted line). The dashed line 
takes only into account the effect of Pauli blocking. 



Next, we turn in Fig. ^ to a comparison of the exact 
calculation for the muon energy spectra (solid lines) with 
the approximation (dashed lines) made on the right hand 
side of Eq. (TJ . Fig. ^ shows the results for the ND and 
Fig. ^ for the FD obtained with the neutrino flux spec- 
ified in Eq. One can see that the approximation is 
in perfect agreement with the exact calculation for muon 
energies i?^ > 1.2 GeV for single pion production and 
Ef_i > 1 GeV for QE scattering. Hence, in this region 
the muon energy spectra are directly proportional to the 
neutrino flux shifted in energy. 

Finally, in Figs. |5^ and|5jD we show the far-near-ratio 
of the muon energy spectra for QE and 1-pion produc- 
tion, respectively. In this case, the dashed lines depict 
the exact result and the solid lines the oscillation prob- 
abilities Pfj_f^ given in Eq. which have been evaluated 
at the shifted energies Eji given in Eq. ©. Again, at 
Ef, > 1 GeV (QE) and > 1.2 GeV (1-pion produc- 
tion) the approximation in Eqs. Q and H1U|) works well 
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FIG. 5: The predicted ratio of tlie far and nearby detector of tine K2K experiment for (a) QE and (b) 1-pion production. The solid 
line is the neutrino oscillation probability. The dashed line denotes the predicted exact ratio of the far and near detector of the K2K 
experiment. 



and a measurement of the far-near-ratio of such pure QE 
or 1-pion event samples in this kinematic region would 
give direct access to the oscillation probability. However, 
as has been discussed above, the observable 1-Ring and 2- 
Ring muon events are superpositions of QE, 1-pion, and 
multi-pion events making the analysis more complicated. 

IV. CONCLUSIONS 

The muon energy spectra of QE and 1-pion events pro- 
vide a complementary approach to experimental extrac- 



tions of the atmospheric muon oscillation parameters at 
K2K. The results are based on quite general kinematic 
considerations and will also be applicable to other fu- 
ture long baseline experiments like J2K, MINOS and the 
CERN-Gran Sasso experiments which plan to use low en- 
ergy beams 0, Q| . Therefore it will be very useful 
to extend this analysis for the beam energy spectra and 
the target nuclei of these experiments. 

Acknowledgment 

J.-Y. Yu wishes to thank the organizers of the ICFP03 
in Seoul for the kind invitation and financial support. 



[1] M. H. Ahn et ai, K2K Collaboration, Phys. Rev. Lett. 

90, 041801 (2003), hep-ex/0212007. 
[2] Y. Fukuda et ai, Super-Kamiokande Collaboration, 

Phys. Rev. Lett. 81, 1562 (1998), hep-ex/9807003. 
[3] C. Y. Walter, Proc. second Intl. Workshop on Neutrino- 
Nucleus Interactions in the Few GeV Region (Nulnt02), 

University of Califonia, Irvine (2002). 
[4] E. A. Paschos, D. P Roy, I. Schienbein, and J. Y. Yu, 

Phys. Lett. B574, 232 (2003), hep-ph/0307223. 
[5] E. A. Paschos, L. Pasquali, and J. Y. Yu, Nucl. Phys. 

B588, 263 (2000), hep-ph/0005255. 
[6] E. A. Paschos and J. Y. Yu, Phys. Rev. D65, 033002 

(2002), hep-ph/0107261. 
[7] E. A. Paschos, Nucl. Phys. Proc. Suppl. 112, 89 (2002), 

hep-ph/0204138. 
[8] P. A. Schreiner and F. Von Hippel, Nucl. Phys. B58, 333 

(1973). 

[9] E. A. Paschos, M. Sakuda, and J. Y. Yu, hep-ph 0308130; 



Published in PRD. 
[10] G. L. Fogli and G. Nardulli, Nucl. Phys. B160, 116 
(1979). 

[11] S. L. Adler, S. Nussinov, and E. A. Paschos, Phys. Rev. 
D9, 2125 (1974). 

[12] I. Schienbein, and J. Y. Yu, Proc. Second Intl. Work- 
shop on Neutrino-Nucleus Interactions in the Few 
GeV Region (Nulnt02), Irvine, Califonia, USA (2002); 
hep-ph/0308010 

[13] Y. Itow, Proc. First Intl. Workshop on Neutrino-Nucleus 
Interactions in the Few GeV Region (NuIntOl), KEK, 
Japan (2001); Nucl. Phys. B (Proc. Suppl.) 112 (2002) 
3. 

[14] P. Lipari, Proc. First Intl. Workshop on Neutrino- 
Nucleus Interactions in the Few GeV Region (NuIntOl), 
KEK, Japan (2001); Nucl. Phys. B (Proc. Suppl.) 112 
(2002) and references there in. 



